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Abstract 

We reinvestigate the conditions for stable matter solutions in the Nambu-Jona- 
Lasinio (NJL) model. In mean field approximation the NJL model can be regarded 
as an extension of the Walecka mean field model to include negative energy fermion 
states. While this extension is necessary to allow for a chiral phase transition, it 
was found some time ago that at the same time it destroys the wanted saturation 
properties of the Walecka model. We reformulate this problem in terms of the 
thermodynamic potential and find that there is indeed a connection between these 
two features. We show that the minimum of the thermodynamic potential which 
corresponds to stable nuclear matter in the Walecka model is shifted from a finite to 
zero effective fermion mass in the chiral NJL model. This shift is closely related to 
the chiral phase transition. Under certain conditions the shifted minima may still 
lead to stable matter solutions but only in the chirally restored phase. We discuss 
a possible interpretation of these solutions as a schematic bag model description. 
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1 Introduction 



One of the major goals in intermediate energy physics is the understanding of strongly 
interacting matter as a function of temperature and density and the transition from the 
hadronic phase to the quark-gluon plasma. Unfortunately nature provides us only with 
few points in the T-p plane, namely the vacuum (T = p = 0) and nuclear matter (T = 0, 
p = p = .17 fm~ 3 ) and much experimental and theoretical effort is necessary to obtain 
information about other regions. A major problem is that in general it is not possible to 
prepare a static system of a given temperature and density. So most points in the T-p 
plane can only be reached for a short instant of time, e.g. on a trajectory of expanding 
matter produced in a heavy ion collision. For zero density and finite temperature the most 
reliable existing "data" therefore come from lattice calculations, indicating that both, the 
chiral phase transition and the confinement-deconfinement phase transition take place at 
the same temperature T c ~ 150 MeV [0. However, because of conceptual difficulties 0, 
there are no firm lattice results at finite densities so far. 

In this situation simplified models may give us at least some ideas about hadronic 
matter at finite temperatures and densities which is necessary for the interpretation of 
the experimental data. Recent results on the lattice |§ and in random matrix models 0] 
indicate that the finite temperature phase transition behaves like a mean-field transition 
at least until a temperature very close to T c . It has been argued that this should still 
be true at finite densities |J. Therefore relativistic mean field models are of particular 
interest. 

Unfortunately until now there is no realistic model of this type which covers the whole 
range from low densities up to the chiral phase transition. The most prominent model in 
the regime near nuclear matter density is the Walecka model || which describes nuclear 
binding and the stability of nuclear matter at saturation density p Q quantitatively as the 
result of the interplay between a large scalar and a large vector field. Many modifications 
of the original Walecka Lagrangian as well as loop corrections have been studied. However, 
because of its simplicity and transparency the linear mean field model (QHD-I) remained 
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the perhaps most important version of the model. 

Although Walecka theory has been originally designed as an effective theory for nuclei 
and nuclear matter at low energies, Walecka-type models are also used for the analysis of 
processes at much higher densities. For instance in ref. |7j an extended Walecka model 
was used to analyze the CERES data at CERN/SPS |§. Here the initial density is 2.5 p a 
at T = 165 MeV. It is clear, however, that at some point the Walecka model must become 
unrealistic as there is no chiral phase transition in the model. In this situation one might 
think of extending the model to negative energy states and to generate the nucleon mass 
dynamically. On mean field (Hartree) level this is equivalent to the Nambu-Jona-Lasinio 
(NJL) model |J with scalar and vector-isoscalar interaction. The chiral phase transition 
in the NJL model has been investigated by many authors, e.g. [0-|L3||. 



It was found, however, by Koch et al. [JL4| that in the NJL model, in contrast to 



the Walecka model, there is no saturation density for stable matter, i.e. the matter de- 
scribed in this model either always expands or collapses, depending on the vector coupling 



strength. (A similar result was obtained by da Providencia et al. [15] who investigated 
the stability of quark droplets in the NJL model within the time- dependent Hartree-Fock 
formalism. They found only expanding droplets.) In order to get the saturation feature 
back the authors of ref. |T^] introduced a new term to the Lagrangian, which effectively 
made the scalar coupling constant density dependent. In section 3 of the present work we 
reinvestigate this problem using the formalism of ref. |ll[] which makes more transparent 
why the saturation properties of the Walecka model get lost in a chiral model. A better 
understanding of this point may help us to construct an improved model in the future We 
will show that the mechanism which stabilizes matter in the Walecka model is basically 
the same as the mechanism which causes a first-order phase transition in the NJL model, 
namely the existence of two minima in the thermodynamic potential. In fact, under cer- 
tain conditions we find stable matter also in the NJL model, however only in the chirally 
restored phase. 

In section 4 we discuss how the NJL phase transition should be interpreted in the 
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light of our results. A general problem of describing the chiral phase transition within 
the NJL model is, that one is dealing with the wrong or at least with an incomplete set 



of degrees of freedom. In ref. |L2] some improvement was made by going beyond the 
standard mean field approach and including mesonic (collective) degrees of freedom. It 
is well known that the low-temperature behavior is dominated by the pion as the lightest 
particle. However, for the finite density phase transition - which is the main task of this 
paper - the more serious question seems to be whether one should deal with nucleons, like 
in Walecka theory or the original NJL model or with quarks like in most NJL papers of 
the past two decades. If the chiral phase transition comes together with the confinement- 
deconfinement one we should have quarks (and gluons) in one phase and hadrons in the 
other. It is obvious that the hadronic phase is not well described by a gas of interacting 
quarks. On the other hand in QCD chiral symmetry is an (approximate) symmetry for 
quarks which does not necessarily translate into a symmetry for nucleons in an effective 
Lagrangian. Ultimately, hadronic matter should be described in terms of hadrons made 



out of quarks, perhaps like in the Guichon model [0, 17] where interacting MIT bags are 
used. In this context our result that stable quark matter in the NJL model does only 
exist in the chirally restored phase, may be seen as an interesting parallel to bag models 
where the hadrons also consist of massless quarks. This will be discussed in more detail. 



2 Formalism 

We consider the following generalized NJL-Lagrangian: 

C = ^(z 7 ^-m )^ + GsK^) 2 + (^7 5 r^) 2 ] - G v {^f . (2.1) 

Here ip is a fermion field with rif = 2 flavors and n c colors. It may be interpreted as 
a nucleon (n c = 1) or as a quark field (n c = 3). Apart from the bare mass m , the 
Lagrangian eq. (2.1) is chirally symmetric (SU(2) L x SU(2) R ). Gs and Gy are positive 
constants of dimension mass -2 , which will be fixed later. 
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Expanding ipip and ip'-y^ip about their thermal expectation values we can derive the 
mean field thermodynamic potential at temperature T and chemical potential /i 1 1 II . 



We restrict ourselves to the Hartree approximation. Furthermore in this paper we only 
consider T = and fi > 0. The result for the thermodynamic potential is 

u MF ^-m^ R ) = w jr> + w jr°(/*) + {m ^f - . ( 2 - 2 ) 



with 



= -(2n f n c ) J ^E p (2.3) 



and 

d 3 p 



(med) 



(fi R ) = -(2n f n c ) J (ja r - E p ) 6{p F - p) , (2.4) 

being the vacuum part and the medium part of the thermodynamic potential (per volume) 
of a free fermion with mass m. E p = \frri 1 + p 2 is its on-shell energy at momentum p. The 



Fermi momentum is given by p F = J fi R — m? B(n R — m 2 ). The vacuum part (eq. (2.3)) 
is strongly divergent and has to be regularized. For simplicity we use a sharp cut-off A in 
three dimensional momentum space. 

In addition to the external parameter /i, ujmf depends on two other parameters, the 
dynamical fermion mass m and the renormalized chemical potential fi R , which are related 
to the scalar density (^V) an d the vector density (ip^ip) at the chemical potential /x by 
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m = m - 2G s (?pip) , 

Hr = fi - 2GV(VV> • 
(2.5) 

These parameters have to be determined selfconsistently by calculating (ipip) an d r) 

from ujmf- It can be shown that the self consistent solutions correspond to the extrema 

of oo mf a s a function of m and fi R . This leads to a set of coupled equations for m and n R 

which reads for T = and \i > 0: 

„ . w f d 3 p m f d 3 p m „. . . . 
rn = m + 2G S (2n f n c )(J - - J ^ - 9( PF - p)) , (2.6) 
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fi R = n - 2G V (2n f n c ) J 9(p F - p) . (2.7) 

Eq. (2.6) is the well-known NJL gap equation for the dynamical mass m at temperature 
T = and chemical potential \i > 0. In general this equation has more than one solution 
and the physical one is the one which minimizes lomf- It is therefore more safe to start 
the analysis from the thermodynamic potential (eq. (2.2)) than from the gap equation 
(eq. (2.6)) because the latter does not tell us which solution is the correct one. However, 
eq. (2.7) has a unique solution for \i R and we can use this equation to eliminate the \ir 
dependence of ujmf- Thus for a given \x we are left with a function which depends on m 
only. More precisely, we define: 

u(p;m) = w MJ? (/i; m, /i B (/i, m)) - w M f(0, 0, m w , 0) . (2.8) 

Here Hr(^, m) is the solution of eq. (2.7) for given /i and m, while m vac is the dynamical 
mass which minimizes the thermodynamic potential of the vacuum. So the second term 
on the r.h.s. is just a constant which shifts the minimum of the vacuum thermodynamic 
potential to zero. Since the absolute magnitude of the thermodynamic potential has no 
physical meaning, we can always subtract such a constant. 

Once we have the thermodynamic potential, other thermodynamic quantities can be 
calculated in the standard way. In particular we have for the baryon number density, 
energy density and pressure 

Pb = — (^V) = ttPf, e = u) + fin c p B , p=-u, (2.9) 

where the energy and the pressure of the physical vacuum was defined to be zero. Then for 
a given \x these formulas have to be evaluated for the values of m and hr which minimize 
the thermodynamic potential. 

For comparison we also consider the QHD-I version of the Walecka model in classical 
mean field approximation. The Lagrangian of this model reads: 

C w = ^[^(idf, - g v uj^) - (m N - g s a)}ip 
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(2.10) 

with a nucleoli field ip, a scalar field a and a vector field u. The classical mean field 
theory which results from solving the equations of motion after replacing the boson fields 
by their space-time independent expectation values is equivalent to the NJL model in 
mean field approximation with n c — 1 and m = (the free nucleon mass), without 
taking into account negative energy states, i.e. neglecting the vacuum part of the free 
thermodynamic potential (eq. (2.3)). Here we made the identification 

1 q 2 1 q 2 

Gs = ^^~ 2 , G v = -^. (2.11) 

So, keeping the NJL notation, the Walecka thermodynamic potential can be written: 

= <rw + {m ~™f - , (2.12) 

and all thermodynamic properties can be derived in the analogous way as discussed above 
for the NJL model. 



3 Stable Matter Solutions 

The saturation of nuclear matter in the Walecka mean field model is usually discussed 
in terms of the energy per nucleon which (with the appropriate choice of the model 
parameters) has a minimum at nuclear matter density p a . In a similar analysis Koch et 



al. |14| found for the NJL mean field that there is no such minimum in the energy per 
nucleon and consequently no stable matter. In order to understand how this difference 
comes about it is useful to recall that the stable nuclear matter in the Walecka model 
is related to a first-order gas-liquid phase-transition with the nuclear matter belonging 
to the liquid phase: If the average density of nucleons is less than p Q the nucleons are 
not uniformly distributed but a mixed phase is formed consisting of nucleon droplets of 
density Pb = Po surrounded by the vacuum with p B = (at T = 0). Now Gibbs' criteria 
tell us that the pressure p* and the chemical potential p* in the droplets must be equal 
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to the pressure and the chemical potential in the vacuum, i.e. p* = and //* < m vac . In 
other words: At some chemical potential fi* < m vac the thermodynamic potential must 
have an additional minimum at m = m* < m vac with 

u)(/i*;m*) = uj(fi*;m vac ) = 0; \m*\ < fi* < m vac . (3.1) 

The existence of this additional minimum is the condition for the existence of stable mat- 
ter. Obviously \i* is just the critical chemical potential for the gas-liquid phase transition. 

In order to illustrate this point, in the left panel of fig. 1 the Walecka thermodynamic 
potential is plotted as a function of the mass parameter m for four different chemical 
potentials. The solid line corresponds to the critical chemical potential /i* where eq. (3.1) 
is satisfied by a minimum at an effective nucleon mass m* = 514.5 MeV. 

As discussed in the previous section the essential difference between the Walecka and 
the NJL mean field is the existence of the vacuum term uJ^ ac ^ (eq. (2.3)) in the latter. It is 
this term which - for a sufficiently large scalar coupling - causes the spontaneous breaking 
of chiral symmetry in the NJL model, i.e. a finite vacuum mass m vac while the bare mass 
m is zero or at least small. At some critical chemical potential the system undergoes a 
phase transition and chiral symmetry gets restored. This chiral phase transition is the 
analogue to the gas-liquid phase transition of the Walecka model: The chirally restored 
phase corresponds to the liquid and the broken phase to the gas. So, according to our 
previous discussion we expect the possibility of stable matter in the NJL model but only 
in the chirally restored phase. 

As an example we consider quarks in the chiral limit. Using standard methods Jl8| we 
can more or less fix the scalar coupling constant Gs and the cut-off A by fitting the pion 
decay constant f n and the quark condensate (uu) in the vacuum. Since the latter is not 
known very accurately we investigate three sets of parameters with different constituent 
quark masses. They are listed in table 1. We do not fix the vector coupling constant Gy 
but treat it as a free parameter and discuss its influence on the results. 
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An important quantity for our later discussion is the bag constant 

B = u(p = 0;m = 0) , (3.2) 

i.e. the pressure difference between the physical vacuum (m = m vac ) and the trivial 
vacuum (m = 0). Its values for the three different parameter sets are also given in 
table 1. 

We start our discussion with parameter set 2 and Gy = 0. The corresponding thermo- 
dynamic potential uj is plotted in the right panel of fig. I for various chemical potentials. 
In the chiral limit Q is symmetric in m. We show only the positive mass region. The 
dotted line corresponds to the vacuum thermodynamic potential (// = 0). The potential 
has a minimum at m — 400 MeV (the vacuum mass we have chosen) and a maximum at 
m — 0. At p — 344.4 MeV the m = maximum turns into a minimum, which becomes 
zero at p* = 378.5 MeV. Finally at p — 409.3 MeV the original minimum disappears. So, 
according to eq. (3.1) we expect stable quark matter at p* = 378.5 MeV with massless 
quarks. This corresponds to a baryon number density of 2.8 times nuclear matter density. 
This is confirmed by the behavior of the energy per baryon number (s/pb) as a function 
of the baryon number density which is shown in fig. 2. The solid line corresponds to 
the massive solutions of the gap equation, the dashed line to the massless ones. The 
energetically favored point is ps = 2.8p G and lies in the chirally restored regime, as we 
expected. 

The phase transition in this example is of first order and takes place at zero density. In 
fact, according to eq. (3.1) it has to be of this type in order to allow for stable matter solu- 
tions. However, the character of the phase transition depends on the model parameters, in 
particular on the vector coupling constant: For parameter set 2 and 0.3Gs < Gy < O.8G5 
there is a first-order phase transition but with a finite transition density (i.e. p* > m vac ), 
for Gy > O.8G5 the phase transition is of second order. In both cases there is no stable 
matter solution. This can be seen from fig. 3 where again the energy per baryon number 
is plotted as a function of density. The left panel corresponds to Gy = \G S) the right 
one to Gy = Gs- Obviously the energetically favored point is ps = in both cases, i.e. 
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the matter is unstable against expansion. 

Our analysis and in particular the comparison between the left and the right panel 
of fig. 1 give some insight why the observation of Koch et al. |14| that the saturation 
properties of the Walecka model get lost when one tries to make it chirally symmetric, 
was not an accidental coincidence. These authors only considered the massive solutions 
of the gap equation and consequently did not find any stable matter in the NJL model. 
The fact, that saturated nuclear matter with massive nucleons does exist in the Walecka 
model, is a direct consequence of its explicitly broken symmetry: The (m — m^) 2 term 
in eq. (2.12) interferes with the m-symmetric w™ 1 *, leading to a minimum somewhere in 
between m = and m = m^. In the chirally symmetric model (NJL in the chiral limit) 
this minimum still exists but is shifted to the symmetric point m = and is identical to 
the minimum which is responsible for the chiral symmetry restoration at finite densities. 

For parameter set 3 we find qualitatively similar results as for set 2: For Gy < O.6G5 
there is a first-order phase transition with zero transition density which leads to stable 
droplets of massless fermions, whereas there is no stable matter for Gy > O.6G5. For 
parameter set 1 the transition density is always greater than zero for positive Gy and 
there is no stable matter solution. 

The main properties of the stable matter solutions are shown in fig. 4 as a function of 
the vector coupling constant Gy. The saturation density ps is plotted in the upper panel, 
while the lower panel shows the binding energy per quark 

E B = -(— m vac ) . (3.3) 

n c ps 

The dashed lines correspond to parameter set 2, the solid lines to set 3. With increasing 
vector coupling the binding energy drops almost linearly. As a result of the increasing 
repulsion the saturation density also drops with Gy, but more slowly. Because of the 
larger bag constant the stable matter is denser and bound more tightly for parameter set 
3 than for parameter set 2. This is also the reason why we did not find any stable matter 
for parameter set 1, where the bag constant was very small. 
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Our results can be understood almost analytically. For m = the energy per baryon 
number is given by 



s(pb) B 3n c ( 3n 



2\ i/3 



Pb Pb A \n f J 
Minimizing this formula we find for the saturation density p$ 

"I (^) 1/3 Pf + G vn y s = B , (3.5) 

provided the energy per baryon number at this point is less than n c m vac , the value at 
Pb = 0. So for Gy = ps behaves like _B 3//4 and it gets reduced by a finite vector 
interaction. The energy density at this point is 

e(p s ) — AB — IGynlpl . (3.6) 

Solving eq. (3.5) for ps until linear order in Gy and inserting the result into eq. (3.6) we 
finally get for the binding energy per quark 

E B ^ rn vac - (—BY* [l + 2Gy(^B)^] , (3.7) 
\n/n c / 37H 

where we neglected terms of order G v . 



For A S> rni ac the following approximate relation holds |19 



Although the cut-off is not very large compared to m vac for our examples, this formula 
works better than 10% for parameter set 1 and even for set 3, where A and m vac are almost 
equal, the error is only 30%. Inserting this formula into eq. (3.7) and taking Gy = we 

can estimate the minimal value of m vac for which Eb > is possible. The result is 

1/2 

5n f n C/ 



m vac Z ( - ) tt/tt ^ 4/ ff , (3.9) 



where we have evaluated the expression for rif = 2 and n c = 3. This result is consistent 
with the fact that we found stable matter solutions for parameter set 2 but not for 
set 1. It is quite remarkable that g = m vac / f n > 4 was also found in ref. POj as the 
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critical condition for solitonic solutions in the NJL model if a sharp 0(3) cut-off is used. 
Reinserting eq. (3.9) into eq. (3.8) we find a minimal bag constant for stable matter 



solutions B > (175 MeV) 4 ~ 125 MeV fm~ 3 , about twice the MIT bag constant |2T 



4 Discussion and Conclusions 

In this section we want to discuss what the results of section 3 mean for the interpretation 
of the NJL phase transition. First let us come back to the question whether nucleon or 
quark degrees of freedom should be used in the model. Starting from low densities, 
nucleons seem to be the more natural choice. However, as outlined above, such a model 
would already fail to describe ordinary nuclear matter. Here the fact that we found a 
stable matter solution in the chirally restored phase does not help us. To the extent finite 
size effects can be neglected, this would mean that nuclei consist of massless nucleons, a 
picture which is certainly unrealistic. If we introduce a small current mass m Q , we could 
obtain a finite effective mass, but still m* would stay well below the Walecka effective 
nucleon mass which is already considered to be too small. So there is no simple extension 
of the Walecka model which keeps its low-density behavior but also allows for a chiral 
phase transition. A possible way out would be to add extra terms to the Lagrangian 
eq. (2.1) as proposed in ref. fl4j] . 



The alternative option is to interpret the NJL model as a quark model, as it is done 
by most present authors. Although it is clear that quarks are not the correct degrees of 
freedom in the hadronic phase, such models agree reasonably well with model independent 
low-density theorems if one uses simple translation rules between nucleon and quark 
quantities (e.g. uiq = mAr/3, a n Q = ct-kn/3 etc.) The hope behind this is that - 

concerning the phase transition - the quarks behave more or less like free, even though 
they are confined in nucleons. On the other hand the stability problem is usually taken less 
serious than in models with nucleon degrees of freedom, since there is no such thing like 
stable quark matter anyway in nature. This point of view is of course very questionable. 
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If we assume that the chiral phase transition for quarks takes place at a finite density 
(i.e. a second-order phase transition or a first-order phase transition with a non-zero 
transition density) our analysis shows that the quarks, if they really behave like quarks 
in an NJL mean field, should repel each other and the matter should tend to expand. 
However, the confinement prevents this expansion, at least for individual quarks. So as 
long as the confinement radius is smaller than the average distance between the quarks, 
i.e. as long as the nucleons do not overlap, the confining forces (or other binding forces 
inside the nucleons) cannot be negligible compared to the scalar and vector mean fields. 

These binding forces are provided by the mean fields themselves if we assume a first- 
order phase transition with zero transition density on the quark level. In this case we 
would find droplets of massless quarks and the density is ps = ps inside these droplets 
and ps = outside. This mixed phase could be interpreted as a schematic picture of a 
bag model description. In a bag model we also have massless quarks at a relatively high 
density inside the bag and zero density outside. 

In fact, in the MIT bag model |%| with massless non-interacting quarks the ground 



state energy of a non-strange baryon is given by [F2TJ 

4vr , 3x-z Q B , , /47r\ 1/3 1/3 , . 

Emit = -j BR + ~ = ~ + ( 3x " z ^ [yJ P * ' (41) 

where R = (j^^) 1 ^ 3 is the bag radius, x = 2.04 and z Q a constant which accounts for 

the zero point motion. So the structure of Emit is exactly the same as that of e/ps for 

massless quarks in the NJL model (eq. (3.4)) if we switch off the vector interaction. This 

is of course not very surprising. Switching off the vector interaction there is no mean 

field left in the chirally restored phase and we are dealing with a gas of free massless 

quarks stabilized only by the bag pressure. So our stable droplets are nothing but large 

quark bags. The difference between eq. (3.4) and eq. (4.1) is that in our thermodynamic 

approach we assumed a very large number of quarks inside a large volume whereas eq. (4.1) 

was derived using the correct boundary conditions for three quarks inside a volume which 

1/3 

can be small. This leads to a different coefficient of the p B term which dominates the 
energy for small radii. However, this difference turns out not to be terribly large. For 
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1/3 

z Q = 2.69 it would even vanish. For z Q = 1.84 as in ref. 21 the coefficients of p B differ by 
about 20%. This corresponds to a 5% difference in the bag radii. As discussed in section 
3 the original MIT bag constant is too small to produce stable matter in the NJL-model. 
The stability densities we find with the larger bag constants of the parameter sets 2 and 
3 (see fig. 4) correspond to bag radii of 0.7-0.8 fm. 

In order to calculate baryon masses the bag energy eq. (4.1) should be corrected for 
spurious center of momentum motion J£| . Again, these correction terms are not present 
in our thermodynamic approach because they are suppressed if the number of particles 
inside the droplet is large. For 3-quark bags, however, they reduce the radii by 20-30%, 
which means the density becomes drastically enhanced. 

Usually a bag model picture of a baryon is considered to be quite different from a 
constituent quark model. However, as just discussed, the NJL mean field may lead to 
the formation of bag-like quark droplets while there are constituent quarks in the model 
as well. Moreover, since the binding energies we found in the previous section are less 
than 15% of m vac (see fig. 4) the bag energies are roughly given by the constituent quark 
relation e/ps — 3m vac . The essential difference to the constituent quark model is that in 
our picture the constituent quark mass is not the mass of the quarks inside a baryon but 
the mass of a single quark state in the vacuum.^ 

These constituent quark states, which do not exist in the bag model are also the reason 
why the NJL model does not confine. For ps — > oo the energy per baryon behaves like 
B/pb for both, the MIT bag as well as for droplets consisting of the massless solutions 
of the NJL gap equation. So in order to increase the MIT bag radius to infinity an 
infinite amount of energy is needed. The same would be true for the NJL droplets if only 



tOn the other hand baryons in the NJL model can also be described by solving Fadeev equations for 
three interacting constituent quarks p4| . In the soliton description ipip strongly varies inside the baryon 
and even changes sign at small radii (see e.g. p5[). So in this model the constituent quark mass is again 
the mass of a single quark in the vacuum rather than in a baryon. It is possible that the question about 
the quark mass inside a baryon is more a question of the most efficient way to sum diagrams (which could 
be different for different processes under consideration) than a meaningful physical question. 
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the trivial solutions existed. However, when the density gets low enough the quarks can 
acquire a constituent mass and only a finite amount of energy is needed to increase the 
radius to infinity. So amusingly the same mechanism which generates the bag pressure, 
namely spontaneous symmetry breaking, prevents it from confining the quarks. In the 
MIT bag model, where the bag pressure is put in by hand, this does not happen. 

Another difference between the NJL model and the bag model shows up when we go 
away from the chiral limit and introduce a small current mass m . Then the bag model 
quarks simply have this mass m . In the NJL droplets, although still much smaller than 
m vac , m* is considerably larger than m . For m = 5.5 MeV this is shown in fig. 5. For 
Gy = m* is about 40 MeV for parameter set 2 and about 20 MeV for parameter set 3 
and it grows for both sets with increasing Gy. Qualitatively this behavior follows from 
the fact that for m ^ the solutions m(ps) of the gap equation do not show a sharp 
transition to m — m but have a long smooth tail. Since the saturation density ps is a 
decreasing function of Gy the corresponding masses m* = m(p s ) increase with Gy. 

In the NJL model the constituent quark mass is closely related to the quark condensate 
in the vacuum (eq. (2.5)) which is the physically more meaningful quantity. Inside the 
baryon droplets the quark condensate is about an order of magnitude lower for m = 5.5 
MeV and even zero in the chiral limit. The baryons are good "vacuum cleaners". If we 
integrate the difference between the quark condensate in the bag and in the vacuum over 
the bag volume we get 

(B\uu + dd\B) = -((#)| M=W - (H>)\p B =o) ■ ( 4 - 2 ) 
Ps 

Here 11 B" stands for "baryon" or "bag". In the chiral limit the r.h.s. reduces to —{ipip)/ps 
which is between 8.0 and 9.3 for parameter set 2 and between 5.3 and 7.7 for parameter 
set 3, depending on Gy. For m = 5.5 MeV we find somewhat smaller values, between 
6.8 and 7.4 for set 2 and between 4.9 and 6.4 for set 3. Although the model is very crude 
these results are in qualitative agreement with the value determined for the nucleon from 
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nN scattering: 

(N\uu + dd\N) = °* N ~ 8, (4.3) 

+ m> d ') 

where we have used \{rn$ + mf 1 ) = m = 5.5 MeV, as before, and the irN sigma term 
a wN = 45 MeV || . 

If the baryon droplets do not further interact, the system will stay in the vacuum- 
droplet mixed phase until the average density reaches the stability density p B = p s . 
Above this density there is no empty space left and the total volume is uniformly occupied 
by massless quarks. This point defines the "macroscopic" finite density phase transition 
compared to the "microscopic" phase transition at zero density which was the starting 
point of this discussion. So the "macroscopic" phase transition would take place at p B ~ 
Ps which can be substantially higher than the critical densities the NJL model predicts 
even for second order phase transitions. For quark degrees of freedom the latter are 
unrealistically low. 

Whereas we could interpret the quark droplets of our thermodynamic approach as 
schematic baryon bags, there is nothing which would correspond to mesonic bags in our 
picture. Mesons have zero baryon number and therefore cannot produce a chirally restored 
bag in their surrounding. However, mesons do of course exist in the NJL model as 
collective quark-antiquark excitations of the vacuum or the medium. In fact, this picture 
is much more powerful than the bag model picture, since it leads to pions as massless 
Goldstone bosons, whereas the MIT bag pion came out too heavy |H] . 



The NJL mesons could mediate interactions between the baryon-droplets. The result- 
ing picture is very similar to the Guichon model fl6], [L7| where MIT bags interact by the 



exchange of scalar and vector mesons. In this way one might obtain a Walecka-like de- 
scription of nuclear matter. Our schematic picture shows a way, how the two ingredients 
of the Guichon model, the bag pressure and the mesonic interaction could be traced back 
to the same origin. 

The main problem to construct such a model quantitatively is the question how to cou- 
ple the mesons in the nontrivial vacuum to the quarks in the bag. The meson properties 
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are quite different inside and outside the bag and one has to take into account boundary 
effects in an appropriate way. This would be the main difference to the Guichon model 
where the mesons do not change their properties at the bag surface. Here a better de- 
scription of the droplet mixed phase which goes beyond the assumption of constant mean 
fields seems to be necessary. 

In summary, we reinvestigated the question of matter stability within the Nambu- 
Jona-Lasinio model in mean field approximation. We confirmed earlier results that the 
NJL model does not allow for stable quark or nuclear matter in the phase with spon- 
taneously broken chiral symmetry but we found stable matter solutions in the chirally 
restored phase. In the model this gives rise to the formation of quark droplets which are 
stabilized by the bag pressure of the surrounding vacuum. These droplets show many sim- 
ilarities to bag model solutions. Although we failed to describe hadronic matter within 
the NJL model we argue that this might be achieved in a second step by considering 
meson exchange between the bags, similar to the Guichon model. 

I would like to thank G.E. Brown for useful discussions. 
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set 1 


set 2 


set 3 


A [MeV] 


650 


600 


570 


G S A 2 


2.14 


2.45 


2.84 


m vac [MeV] 


313 


400 


500 


U [MeV] 


92.6 


93.9 


92.5 


(mm) 1 / 3 [MeV] 


-248.9 


-244.9 


-242.6 


B [MeV/fm 3 ] 


73.9 


135.3 


224.8 



Table 1: Three sets of model parameters (0(3) cut-off A and scalar coupling constant 
Gs) and the corresponding vacuum properties for the NJL model with quark degrees of 
freedom in the chiral limit: constituent quark mass m vac , the pion decay constant the 
quark condensate (uu) and the bag constant B. 
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Figure Captions 



Fig. 1 Thermodynamic potential w as a function of the mass parameter m for different 
chemical potentials. 

Left panel: Walecka model with m N = 939 MeV, m s = 550 MeV, m v = 783 MeV, 
g s = 10.3 and g v = 12.7. The different lines correspond to p, — (dotted line), 
p = 915.1 MeV (dashed line), p = p* = 923.15 MeV (solid line) and p = 939.65 
MeV (dashed-dotted line). 

Right panel: NJL for parameter set 2 (see table 1) and Gy = 0. Here the different 
lines correspond to p — (dotted line), p = 350 MeV (dashed line), p — p* — 378.5 
MeV (solid line) and p = 410 MeV (dashed-dotted line). 

Fig. 2 Energy per baryon number e/pg as a function of the density ps for parameter set 2 
and Gy = 0. The solid line corresponds to the massive solutions of the gap equation, 
the dotted line to the massless ones. 

Fig. 3 The same as fig. 2 but for Gy = \G S (left panel) and Gy = G s (right panel). 

Fig. 4 Saturation density ps (upper panel) and binding energy Eb (lower panel) of stable 
quark matter as a function of Gy/Gs- The dashed lines correspond to parameter 
set 2, the solid lines to set 3. There is no stable quark matter for parameter set 1 
in the positive Gy regime. 

Fig. 5 Effective quark masses m* for stable quark matter as a function of Gy/Gs- The 
dashed (solid) line corresponds to the scalar coupling constant Gs and the cut-off 
A of parameter set 2 (3), but using a finite current quark mass m = 5.5 MeV. 



